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ABSTRACT 

Aims. Molecular clouds act as primary targets for cosmic-ray interactions and are expected to shine in 7-rays as a 
by-product of these interactions. Indeed several detected 7-ray sources both in HE and VHE 7-rays (HE: 100 MeV < 
E < 100 GeV; VHE: E > 100 GeV) have been directly or indirectly associated with molecular clouds. Information on 
the local diffusion coefficient and the cosmic-ray population can be inferred from the observed 7-ray signals. In this 
work we explore the capability of the forthcoming Cherenkov Telescope Array Observatory (CTA) to provide such 
measurements. 

Methods. We investigate the expected emission from clouds hosting an accelerator, surveying the parameter space for 
different modes of acceleration, age of the source, cloud density profile, and cosmic-ray diffusion coefficient. 
Results. We present some of the most interesting cases for CTA regarding this science topic. The simulated 7-ray 
fluxes depend strongly on the input parameters. In several cases, we find that it will be possible to constrain both the 
properties of the accelerator and the propagation mode of cosmic rays in the cloud from CTA data alone. 

Key words, astroparticle physics - radiation mechanism: non-thermal - ISM: clouds - cosmic-rays - gamma rays: ISM 



1. Introduction 

Emission in HE- VHE 7-rays is expected in spatial coinci- 
dence with molecular clouds, resulting from the hadronic 
interaction between cosmic-ray (CR) particles and the 
dense material in the cloud acting as a target. Indeed, 



GeV and TeV domain (see, e.g., 


Aharonian et all 


2008 


Albert 2007: Aharonian et all 2008 


; Giuliani et al. 


2011 


Ackermann et al.ll2012t lAharonian 2012). Moreover, it has 



been suggested that some of the as yet unidentified 7-ray 
sources might also be MCs illuminated by CRs that es- 
caped from a n accelerator located inside the cloud or in 
its proximity (iMontmerld 19791: Aharonian 



Atovan|[T996t 



iGabici et aLll2007l : iRodrfguez Marrero et al.ll2008l ). In such 
cases the modeling of the emission involves the parametriza- 
tion of the diffusion of charged particles. The diffusion co- 
efficient is in general conside red to be energy-dependent 
(|Ginzburg fc Svrovatskii|[T964D , with lower energy particles 
diffusing more slowly than higher energy ones, under the 
same medium conditions. When those diffused charged par- 
ticles (protons or heavier nuclei) interact with target ma- 
terial in a density enhancement of the interstellar medium, 
such as a molecular cloud located in the vicinity of the accel- 
erator, significant 7-ray emission is expected due to the pro- 
duction and subsequent decay of neutral pions. 7-ray emis- 
sion produced in massive molecular c l ouds was predicted 
long ago (see e.g . iBlack fc Faziol 119731: iMorfAll et all 119841: 
[Ah aroma 3 ll99lf ). The study of this emission is extremely 



useful in unveiling the physics of CR sources. Due to energy 
dependent propagation effects the 7-ray spectrum from the 
molecular clouds may differ significantly from the spectrum 
observed at (or closer to) the accelerator. This may explain 
discrepancies in the particle spectral indeces inferred from 
the same source at different frequencies, even if all particles, 
leptons and hadrons, are accelerated to the same power-law 
at the source. In this scenario, a great variety of 7-ray spec- 
tra is expected, depending on several parameters, includ- 
ing: the age of the acceleration, the distance between the 
cloud and the accelerator, the duration of the injection of 
CRs and their diffusion coefficient. This can produce a vari- 
ety of different GeV-TeV connections, some of which could 
expla i n the observed p henomenology (see, e.g.. iFunk et al.l 
120081: iTam et aLlfeOlOD . Injection and propagation of CR 
have been recently put forward to explain a number of the 
TeV sources currently known, especially those for which 
there is no, or there is a s patially displaced, GeV coun- 
terpa rt fe.g. seelFuiitall2009l;lLi fc 'Chenll2010t lG abici et al.l 
l2010l : lTorres et al.ll2010l : lOhira et al.ll2011h . 



Given the expected CTA angular resolution and sen- 
sitivity, variations in flux at less than 5 pc bins at 5kpc 
distances could be resolved. Changes in the CR spectrum 
could be derived accordingly, leading even to, e.g., the 
derivation of a diffusion coefficient as a function of energy 
and position (D(E,r)). The measurement of such spatial 
variability in the diffusion coefficient would be an impor- 
tant result in CR physics. 
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In this paper, the 7-ray emission due to an accelerator 
inside a molecular cloud is calculated. The expected CTA 
measurement of such emission is then derived taking into 
account the simulated CTA response functions. The CTA 
Observatory is described in Section 2. The calculation of 7- 
ray emission and the physical parameters of the scenario are 
described in Section 3. The simplified case of flat density of 
the target material (i.e. flat density profile of the molecular 
cloud) is investigated in Section 4 along with the CTA ca- 
pabilities in distinguishing the parameter space. Small and 
nearby clouds are investigated in Section 5. Section 6 deals 
qualitatively with a more realistic case of a peaked density 
profile. Conclusions are given in Section 7. 

2. The Cherenkov Telescope Array (CTA) 

CTA is an international project for the development of 
the next g eneration ground-based 7-ray instrument (see 
lActis et alj|2011l) . The detection of 7-rays (E >10 GeV) 
with ground-based facilities is possible thanks to the imag- 
ing atmospheric Cherenkov technique (|Weekesll2003[ ) . VHE 
7-rays interact with nuclei in the atmosphere producing 
a cascade of particles, where velocities are larger than 
the speed of light in the medium, leading to Cherenkov 
light emission. The resulting Cherenkov light flashes may 
be imaged by Imaging Atmospheric C herenkov Telescopes 
(IAC Ts) (for a recent review, see iHinton fc Hofmaiml 
I2009T) . The shower images are then used to reconstruct the 
energy and direction of the original particle. Particle cas- 
cades can also be initiated by CRs and constitute the main 
source of background. In this case the showers are generally 
broader than those initiated by primary photons. Gamma- 
hadron separation can be achieved thanks to the differences 
in size and shape of the shower images. The showers illu- 
minate a pool at the ground level. The radius of the light 
pool depends on the height of the observatory and on the 
energy of the primary photon. An array of IACTs allows 
for better sampling of the Cherenkov light distribution of 
a given event. From a stereoscopic view of the same event, 
the reconstruction of the direction of the primary photon 
and the background rejection are improved with respect to 
a stand-alone telescope observation. 

CTA will significantly advance on the present generation 
IACTs: it will feature an order of magnitude improvement 
in sensitivity at the core energy range of 1 TeV, improve 
in its angular and energy resolution, and provide wider 
energy coverage, see lActis et al.l (1201 ll ). Indeed, the array 
is expected to have an unprecedented sensitivity down to 
~50 GeV and above ~50 TeV, establishing a strong link 
to the satellite-based operations at low energies, namely 
the Large Area Tele scope on board the Fermi satellite, sec 
I At wood et alJ (|2009T) and water Che renkov experime nts at 
the highest energies (c.g, HAWC, see lGoodmar3l201Cfl . The 
gain in sensitivity is due to the increase in the number of 
telescopes. The widening of the explored energy range is 
due to a combination of different-sized telescopes in dif- 
ferent parts of the light pool. Large size telescopes (LST, 
with a dish of ^23 m) will be placed at the center of the 
array. Thanks to their large mirror area, dim flashes from 
the low energy events (~50 GeV) are expected to be recon- 
structed. Tens of medium size telescopes (MST, with a dish 
of ~11 m) will be placed in a surrounding ring, covering a 
large fraction of the light pool and thus enhancing the re- 
construction of medium energy (^1 TeV) events. Finally, 



the outer regions will be composed of small size telescopes 
(SST, with a dish of ^7m) enlarging the effective area of 
the array for the bright but rare high energy events (above 
~50 TeV). Both a southern and a northern hemisphere ob- 
servatory are foreseen. 

3. An accelerator inside a cloud 

If a power-law energy spectrum (J p (E p ) = RE" 1 ) is as- 
sumed for the intensity of primary CRs, the resulting 7- 
ray spectrum due to hadronic interactions would also fol- 
low a power-law spectrum (F(E) oc E~ r ). However, if we 
consider an energy-dependent diffusion coefficient, the CR 
spectrum may differ from a simple power-law near the ac- 
celeration site. The spectrum of the accelerated CR can 
be expressed as J p (E p ,r 7 t) = (c/4ir)f, where f{E pi r 1 t) is 
the distribution function of protons at time t and distance 
r from the source. The dis tribution function satisfies the 
diffusion-loss equation (e.g.. iGinzburg fe Svrovatsk"iHll964l ) 
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where P = —dE p /dt is the continuous energy loss rate 
of the particles, Q = Q(E p ,r,t)5(R) is the source func- 
tion (for injection), and D(E p ) is the diffusion coefficient. 
Here, we assumed that the source is point-like and lo- 
cated at the origin of the coordinate system. Solutions 
to this equation have been extensively studied for differ- 
ent cases, considering either spatially constant diffusion 
coefficient (lAtovan et al.lfl995t lAharonian fc Atovanl 119961 : 



iRodn'guez Marrero et all 120081: [Gabici et al.1 I2009D or no 

CR accelerator near the cloud, i.e. Q=0, (passive clouds 
where the only 7-ray emi ssion arises from t he contribu- 
tion of the CR background, iGabici et al.ll2007T) . We investi- 
gate the case of an accelerator positioned at the center of a 
molecular cloud. This is an idealized case, but it allows to 
study the impact of an enhancement of CR content, above 
and beyond the passive cloud case. The 7-ray emission can 
be calculated in concentric shells of increasing radius, each 
shell retaining the footprint of the diffusion coefficient and 
of the cloud density. The study of such footprint is done for 
a simple symmetrical and homogeneous system, where ex- 
pected spectral and morphological behaviors can be shown 
clearly 

The acceleration and diffusion processes are com puted 
following the approach of lAharonian fc Atovanl (j!996l) . The 
diffusion coefficient is assumed to depend on the CR en- 
ergy only, as: D(E p ) = D 10 {E p / 10 GeV) 5 cm 2 s _1 . More 
details of the flux calculation are given in the Appendix [5] 
The resulting flux in 7-rays is mainly dependent on the dif- 
fusion coefficient, the age of the accelerator, the type and 
spectrum of injection of accelerated particles, and on the 
density and mass of the cloud. An impulsive source of par- 
ticles corresponds to the case when the bulk of rclativistic 
cosmic-rays are released during times much smaller than 
the age of the accelerator itself. When the timcscales are 
comparable, the source is referred to as a continuous injec- 
tor. All the parameters are free and might assume slightly 
different values to those studied here. The intervals for the 
values of the parameters adopted in this work are given 
below: 



— Diffusion coefficient 
[10 26 ..10 



281 cm 2 s _1 ); 



Slow to fast (e. 
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— Diffusion coefficient energy dependence: 6 = [0.3..0.6]; 

— Age of the accelerator: [10 3 .. 10 5 ] years 

— Type of accelerator: Impulsive / Continuous 

— Spectrum of injection: 7= [2.0. .2. 5]; 

— Fraction of energy in input (total energy in form of CRs 
W p = 77IO 50 erg): r\ = [0.3. .3]. 

However, we mainly concentrate, as an example, on the 
case were the injection slope is 7 = 2.2 and the diffusion 
coefficient dependence on energy is 5 = 0.5. This pair of pa- 
rameters satisfies the observed CR spectrum data. Indeed, 
with these values, the index of the equilibrium spectrum in 
the galaxy is expected to be 7 + S — 2.7. The typical value 
of th e diffusion coeffic ient in the galaxy is Diq = 10 28 cm 2 
s _1 (jBerezinskv 1990). However, this value is very uncer- 
tain and depends on the level of the magnetic turbulence 
in which particles propagate. For example, higher level of 
turbulence will in general lead to a suppression of diffusion. 
The total energy input is taken as W p = 10 50 erg (77 = 1) 
in the impulsive case, while the energy injection rate in 
the continuous case is of L p = 10 37 erg s _1 , resulting in 
the same total input for accelerators with an age of a few 
hundreds of thousand years. 

The spectra of 7-rays produced b y proton-proton inter- 
action s have been computed following lAharonian &: Atovanl 
(|l996l ). where a delta- function approximation has been used 
to model the interaction cross section. 



3.1. CR sea penetration in the cloud 

The CR background (also referred to as 'sea') is ubiquitous 
in the Galactic plane. For the energy range considered here, 
it is assumed to be well described by the locally meas ured 
differential spectrum as (e.g., see lAdriani et al.ll20lih 



-2.7 



J e (E p ) ~ l.5E p - GcV cr. 



- 2 -- 1 sr^ 1 GeV- 1 . 



(2) 



Therefore, to compute the total CR content inside the 
cloud, it is necessary to calculate the degree of penetration 
of the CR sea in the simulated cloud. The comparison of 
relevant timcscalcs (pp lo sses and diffusion t imescale) has 
been done following, e.g.. iGabici et al.l (|2007t l. The energy 
loss timescale for proton interactions is 



1 



'pp 



(3) 



where «h is the density of the gas, c is the speed of light, k 
is the inelasticity (assumed to be k ^0.45 throughout the 
paper) and a pp = 33 mb is the cross-section of the pro- 
cess. a pp is mildly dependent on the energy of the particles, 
however, for the energies considered h ere, the assumption of 
energ y independence is satisfactory (jAharonian fc Atovanl 
|1996|) . The diffusion timescale, i.e. the time it takes a CR 
to diffuse from the edge to the centre of the cloud, can be 
expressed as 



R: 



T diff 



loud 



6D(E P )' 



(4) 



where i? c ioud = 20pc is the radius of the cloud in the ex- 
ample investigated here. For a mass of 10 5 M Q , this results 
in a uniform density of njj = 130 cm -3 , rather typical for 
a molecular cloud. The timescale of Eq. 2J in the impulsive 
case, represents the time at which the maximum of particle 



flux i s reached at the distance i? c ioud (jAharonian fc Atovanl 
Il996l ) . Penetration is ensured for particles of energies higher 
than the energy resulting from equating Eq. ([3]) and Eq. (U). 
For the smaller value of the diffusion coefficient adopted 
here (i.e. Dio = 10 26 cm 2 s _1 ), the minimum energy of 
CR that can penetrate fully the cloud is E p w 2 GeV. 
This is shown in Fig. [1] as the intersection of the dashed 
line representing the pp loss timescale and the black solid 
line representing the diffusion timescale in the case of slow 
diffusion. For faster diffusion, the minimum energy is even 
lower. This conclusion does not depend on the preperties of 
the accelerator, so as to say, it is valid both for continuous 
and impulsive accelerators, but depends on the diffusion 
coefficient, the density of the gas and the size of the cloud. 
Therefore, the CR background density is always added to 
the CR spectrum from the accelerator itself in the cases 
shown in Fig. [T] From the latter figure, thus, it is possi- 
ble to derive the degree of penetration to the core of the 
cloud for CRs of different energies, assuming a flat density 
distribution. 
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Fig. 1. Comparison of relevant timescales calculated via 
Eq. (|3]) and Eq. (01 . The dashed line represents the loss 
timescale through pp channel. The shaded areas represent 
the diffusion timescale for the interval 8 — [0.3..0.6]. The 
solid lines are for S = 0.5. The crossing of the dashed and 
solid lines represent the minimal energy at which total cloud 
penetration is fulfilled for a diffusion coefficient of Z?io = 
10 26 cm 2 s" 1 (black), D w = 10 27 cm 2 s" 1 (red), D w = 10 28 
cm 2 s _1 (green). 



4. CTA response 

4.1. Detectability 

The 7-ray emission has been simulated for a molecular 
cloud with a mass of M 5 = 10 5 M Q and a radius of 20 pc 
(hence with an average density of nn = 130 cm~ 3 ) located 
at a distance of d =1 kpc, with the accelerator parameters 
given above for a continuous and impulsive source. Giant 
molecular clouds as close as ~1 kpc distance might be un- 
common, however for a list see iDame et al.1 (|1987h . The 
angular extension is ~ 1 deg in radius. Most of the results 
can be rescaled to an arbitrary cloud mass and distance by 
recalling that, for a given CR intensity in the cloud, the 
expected gamma ray flux is expected to scale as oc A/5/d 2 
and the cloud apparent size as oc R c \ ud/d. 7-ray fluxes were 
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calculated for a permutation of age, acceleration mode, and 
diffusion coefficient as outlined in Section [3] Fig. [2] shows 
that for such permutations, all the cases can be detected 
with 50 hrs of CTA observations except for the case of an 
old impulsive accelerator and fast diffusion. The sensitivity 
of the instrument is scaled for the extension of the source, 
as detailed in Sec. 14.21 The calculated integral fluxes are 
very similar for some of the permutations. For example, in 
the case of fast diffusion and continuous acceleration, the 
calculated flux is constant with age (red, green and blue 
squares at D 10 = f0 28 cm 2 s _1 ). This can be explained 
by the fact that the bulk of the particles contributing to 
7-ray emission at energies above 10 GeV can diffuse out of 
the cloud in a time smaller than the age of the accelerator. 
Therefore only the latest generation of accelerated parti- 
cles contributes to the 7-ray signal and a steady state is 
reached. 
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Fig. 2. Energy flux, for the set of permutations described in 
the text. Filled triangles are for impulsive cases and squares 
are for continuous acceleration (accelerator age of 10 3,4,5 
years are represented by red, green, blue markers, respec- 
tively). The dotted line re presents the CTA integral sensi- 
tivity above 50 GeV as in lActis et all (|2011h . scaled for 1 
degree radial extension. 



Once a spectrum and source morphology are obtained 
using CTA these would act as diagnostic tools with which 
to reconstruct the initial parameters of the accelerator. 

4.2. Spectral features 

A constraint on the diffusion coefficient can come from the 
identification of a break in the 7-ray spectrum integrated 
from the entire cloud region. The break can be related to 
the minimum energy that can diffuse in the entire cloud 
over a timescale comparable to the age of the accelerator. 
By equating Eq.|4j which gives the diffusion timescale, with 
the age of the accelerator, one obtains the constraints shown 
in Fig. [3J which correspond to: 



E, 



p, break 



10 



loud 



1/5 



GeV. 



(5) 



= 10 28 cm 2 s _1 ) and 



For scenarios with fast diffusion (-D10 

energy dependence parameter in the range S = [0.3. .0.6], 
the corresponding break in 7-ray emission will always be at 
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Fig. 3. Boundaries for a break at energies higher than 70 
GeV (in the emitted 7-ray spectrum) for different ages 
of the accelerator (10 3,4 ' 5 years are represented by solid, 
dashed, dotted curves, respectively). The parameter space 
at the left of the boundaries results in a break in the spec- 
trum within the energy domain of CTA. 
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Fig. 4. CTA expected performances in the reconstruction 
of the intrinsic spectral model. Top panel: £>io = 10 26 cm 2 
s _1 , 6 = 0.4, 7 = 2.3, J) = 1/3. Bottom panel: D 10 = 10 28 
cm 2 s _1 , S = 0.6, 7 = 2.1, r] = 1. The blue points illustrate 
one of the possible realizations of the simulated spectral 
points from 50 hours of CTA observation time, the red line 
identifies the intrinsic spectrum and the black lines repre- 
sent the accepted models (x 2 /dof +1). The red intrinsic 
spectrum is superimposed to one of the accepted models, 
which, for the case shown in the top panel, is superimposed 
to the only accepted model. 



energies below the CTA energy acceptance. This is accurate 
for impulsive accelerators. In the continuous acceleration 
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case, new injections of high energy particles will smooth 
the effect of a break. At energies higher than the break 
given by Eq. [5J the particle spectrum will follow a power- 
law form composed of the slope of the injection spectrum 
and the energy dependence of the diffusion coefficient (i.c 
7 + S). Therefore the 7-ray emission will show a power- 
law behavior, thereby reducing the ability to constrain the 
parameter space from 7-ray data in the cases when -Ebreak 
is below 70 GeV. 

The reconstruction of a break is indeed a powerful tool. 
If a break is not present, and the CTA detected spectrum is 
a simple power- law, a plethora of models will fit the spectral 
points. An example is given in the following, where the 
parameters are searched for in a grid for the intervals given 
in the Appendix [A"l 

Let us assume that a molecular cloud with measured 
mass and distance is detected at TeV energies. Let us fur- 
ther assume that from multiwavelength observations we 
identified a possible accelerator of CRs responsible for the 
gamma ray emission and that an estimate of the age of 
that accelerator is known. In Fig.@]wc show the simulated 
spectra for the gamma ray emission for such an acceler- 
ator, which is assumed to be impulsive and with an age 
of 10 4 years. The other parameters are varied on a grid 
(see Appendix |X| and the corresponding observed spec- 
trum is simulated from CTA "responses". By using a x 2 
minimization procedure, the best fitting model belonging 
to the parameter space grid is found, along with a repre- 
sentative sample of the models within the x 2 /cfc>/ + 1 con- 
tour (do/=degrees of freedom), which is those models with 
a x 2 /dof not in excess of unity from the best fit model. The 
top panel of Fig. [4] shows the case of D w = 10 26 cm 2 s _1 , 
S = 0.4, 7 = 2.3, T] = 1/3. Thanks to high flux reached in 
this case and the presence of a break, the spectrum is recon- 
structed easily to the intrinsic parameters, with a break in 
the 7-ray spectrum at E ps 2 TeV and slopes Ti = 2.3 and 
T2 = 2.7 below and above the break, respectively However, 
for lower fluxes, the reconstruction can be more uncertain. 
This is shown in Fig. U (bottom panel), where it can be 
seen that the intrinsic model (-D10 = 10 28 cm 2 s _1 , S = 0.6, 
7 = 2.1,J7 = 1) does not even provide the best fit. 

Additional and more stringent constraints can come 
from studying the spatial dependence of the 7-ray spectrum 
from the inner to the outer region, which depends on the 
diffusion inside the cloud. We slice the expected emission, 
projected onto the sky, in concentric shells (see Appendix). 
We choose the linear size of the shells to match the expected 
angular resolution of CTA at the lowest energies resolvabl e 
by the array (e.g. -0.25° at 50 GeV, see lActis et aDl201lh . 
with 5 concentric shells, with shell 1 being the closest to 
the accelerator. Fig. [S] shows an example of the predicted 
fluxes and Fig. [5] shows the expected CTA energy spectra. 

Fig. [5] shows also the expected CTA sensitivity scaled 
with the size of the source. In order to obtain the CTA sen- 
sitivity for an extended source, t he point-source sensitivity 
is taken from lActis et al.l (|201lf ) and then scaled with an 
appropriate energy dependent factor. This factor is related 
to the optimal cut on the angular size of the source. The ra- 
tio between sensitivity for a point-like source (PS) and an 
extended (EXT) source is ©ext/0ps- In Fig. [5] we show 
the PS sensitivity for CTA together with the scaled one for 
©ext- We show this to illustrate the capability of CTA. 
However, to calculate the spectral points and profiles from 
CTA simulated observations, we followed the procedure de- 



tailed in IBernlohr et al.l j(2U 1 2 : . To evaluate the spectral 
CTA response for each shell, we simulate the background 
counts coming from a region as extended as the outer bor- 
der of each shell. Contamination from adjacent shell is not 
taken into account. 

The LAT instrument on board of the Fermi satellite 
will provide at least 5 years worth of data by the time 
that the full CTA array will be in operation. Figs. [5] and 
[6] show the expected 5 year point source sensitivity for 
Fer mi /LAT. This is ca lculated from the 1 year sensitivity 
in (jAtwood et al.l[2~009f ). linearly scaled with time at high 
energies (> 10 GeV). The linear scaling is expected in the 
signal limited regime. 

The simulations show that the total cloud emission, as 
well as that from some individual shells, are well above the 
limit of detection with CTA. For the particular case in Fig. 
IH1 the low energy particles still have not diffused out to the 
outer shells (the last shell is in magenta in the plot). Indeed 
the concave shape of the spectrum evidences the dominant 
contribution from the CR background at low energies on 
the outer shells. The concavity in the spectrum is expected 
at energies below the energy range of CTA. However, it will 
be possible to distinguish a hardening of the spectrum for 
shells with increasing distances from the accelerator. These 
features can be seen for middle age accelerators (age=10 4 
years) only for Diq = 10 26 cm 2 s . For faster acceleration, 
the spectra of all shells will be a simple power-law with 
same photon index, thus not distinguishable. 

4.3. Morphology 

The different scenarios can also be disentangled by inves- 
tigating the morphology and extension of the emission re- 
gion. We used the fluxes and shell area from the scenarios 
described above in order to simulate an excess map from 
the simulated r espon se functions of CTA, as detailed in 
IBernlohr et al.l (|2012t) . We then create a profile from the 
excess map, integrating azimuthally the counts for bins in 
increasing angular distance from the center of the map. The 
profiles are then weighted by the bin area. The shape of such 
profile depends on the parameters of the simulated scenario, 
an example of which is shown in Fig. [71 These profiles are 
easily distinguishable from one another. Extensions of the 
7-ray emission depends on the parameters studied here, 
with some general trends. Older sources are always more 
extended than younger sources, as the lower energy parti- 
cles will have diffused further from the center of the cloud 
and thus from the accelerator. Emission due to continuous 
accelerators will present steeper profiles due to the freshly 
accelerated particles in the center of the source. Faster dif- 
fusion also leads to a larger extension. The profiles in Fig. 
[Jjarc normalized to their respective maximum counts, with 
flux decreasing with age in the impulsive case and oppo- 
site behavior in the continuous case. However we do not 
expect many of these high flux objects, therefore it is use- 
ful to predict the maximum distance at which an object 
is expected to be detected and resolved, depending on its 
intrinsic luminosity. 

The maximum distance (<i ma x) at which a source can 
be detected depends on the sensitivity of the instrument 
and on the angular size of the source, 8 SIC = £,R c ioud/d (C 
represents the fraction of the cloud radius that contains 
the 68% of the counts in profiles as shown in Fig. [7]) . For 
a source of constant luminosity, e? 2 lax = L SIC / (4irF seas ), 
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Fig. 5. The 7-ray flux of the 5 shells considered. The color code goes from the inner to the outer shell: red, green, 
blue, yellow, magenta. Th e black line is the total flux from the molecular cloud, i.e. the sum of the 5 shells, and can be 
compared with Fig. 10 of lAharonian &: Atovanl (|1996l ). The two panels show the example for a continuous accelerator, 
with Dio = 10 26 cm 2 s _1 (slow diffusion, left) and Dio = 10 28 cm 2 s _1 (fast diffusion, right). In the case of fast diffusion, 
the spectrum is similar for all the shells, in particular the blue and yellow lines are coincident. The accelerator age is 10 4 
years. 




Energy E(TeV) Energy E(TeV) 

Fig. 6. CTA response to the scenario with Dio = 10 26 cm 2 s _1 and an age of the accelerator of 10 4 years; both for 
continuous acceleration (left) and impulsive acceleration (right), investigated for a 50 hours integration time. The color 
code follows that of Fig. [3J 
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Fig. 7. Profiles of the photon count, normalized to the re- 
spective maximum. Here we show the example of impulsive 
(Top) and continuous (Bottom) acceleration, Dio = 10 26 
cm 2 s _1 , at an accelerator age of, from bottom to top, 10 3 
(red), 10 4 years (green), and 10 5 years (blue). Error bars 
are set to 10% of the count number, to mimic the expected 
error on the effective area for the array used in the counts 
determination. 



where L s 



47rF src G? 2 . c is the isotropic luminosity of the 



source. The sensitivity of the array, for an extended source, 



is F scns (d) = F PS 8 sl . c /6ps = F PS /8 PS x £,R c i oud /d, where 
Fps is the expected point-source sensitivity of the array 
and 8ps its angular resolution. For CTA, the point-source 
sensitivity is F ps ,cta(> 50GeV) « 10- 12 TeV/cm 2 /s and 
Aps.ota ~ 0-25° its angular resolution, at low energies 
(jActis et al.ll2011f ). Therefore, for F sens (d max ), we will have: 



14 



4 x 10 33 ergs _1 
lOpc 



F] 



PS, CTA 



PS 



cps 
0.25 c 



loud 



kpc, 



(0) 



The maximum distance at which a source can be re- 
solved is when the angular extension of the cloud equals 
the PSF of the instrument (# src = 8ps), therefore the max- 
imum distance is 



-<<rcs,max 



G 



D 



4 x 10 33 ergs~ 



1/2 



/ Fps,CTA 



PS 



1/2 



kpc. (7) 



At this distance we will be able to discriminate between 
the point-source and extended case, but the exact deter- 
mination of the extension might need longer observation 
times. Fig.[5]shows the evolution of the maximum distances 



6 



G. Pedaletti et al: CTA and cosmic ray diffusion in molecular clouds 



of detection and rcsolvability as a function of the luminosity 
of the source. As an example, let us consider a passive cloud, 
i.e. a case where only the CR background is included. With 
this assump tion, the expected signal is (|Aharonianl 1 1 9 9 ll : 
IGabicil[200l : 



/ 77i 

F(> B 7 ) ~ l x 1(T 13 k ' 



-1.7 



VTeV, 

/ M 



Vio 5 M Q y vikpc 



— 9 — 1 

cm \ 



(8) 



where k is the enhancement factor of CRs, assumed to 
be unity for passive clouds. Therefore a cloud of 10 5 M© 
would be detected out to only ~ 1 kpc, due to its expected 
isotropic luminosity of i; so k 2x 10 32 erg s _1 . If such cloud 
was to be more compact than the ones investigated here, 
the horizon of its detectability would be larger, according 
to the scaling given in Eq. [S] It has to be noted that enhanc- 
ing the CR content of the cloud by a factor n > 20 would 
allow the detection of a cloud of such mass in the entire 
galaxy. For comparison, k ~ 20 would correspond to one 
of the cases exemplified in Fig. [21 specifically a continuous 
accelerator of 10 4 year age and Z?io = 10 27 cm 2 s _1 . 



o 

D. 

M 

o 
a 

C3 



X 




horizon of detectability 
horizon of resolvability 



10" 



valid only for an infinite field of view with flat acceptance, 
while wc expect a degradation of sensitivity with increasing 
angular dista nce from the center of the field of view (for a 
discussion see iDubus et al.l l2012t) . Many of the clouds are 
very extended with respect to the expected angular accep- 
tance of CTA and would need either to be scanned or to be 
observed in divergent mode, therefore increasing the obser- 
vation time for a mapping of the entire object (fo r a general 
survey in divergent mode see lDubus et al.|[2012ft . Expected 
fluxes and the corresponding scaled sensitivity are given in 
Fig.[F3 The ranges in extension refer to the cloud boundaries 
in latitude and longitude. 




Perseus OB2 Rho Oph 



Cepheu^ 



OrionB 
OrionA 



10 



12 14 

[degl 



Fig. 9. Expected flux from the GMC in the Gould Belt vs 
their extension (e xtension is 70% of the boundaries listed in 
iDame et al]|1987l ). Expected flux (> 100 GeV) is calculated 
through Eq. [51 with k = 1. The solid line is the on-axis sen - 
sitivity of CTA (> 100 GeV, 50 hours; lActis et alJl201lh . 
linearly scaled with the extension of the source. This opti- 
mistic scaling holds for an infinite field of view with flat ac- 
ceptance, while we expect a degradation of sensitivity with 
increasing angular distance from the center of the field of 



Luminosity (erg s ) 



Fig. 8. Horizons of detectability (solid) and rcsolvability 
(dotted) as a function of isotropic luminosity (£ = 0.5, other 
parameters as in text). 



5. Passive clouds: giants and cloudlets 

Giant molecular clouds (GM C) as close as ~1 k pc distance 
are uncommon (for a list see IDame et al.lll987h . These arc 
passive clouds (i.e. only the CR background is included) 
that would also be interesting for detection. The feasibility 
of a detection depends mainly on their mass, distance and 
most of all e xtension of the gam ma-ray emission, as already 
discussed in lAcero et al.l (|2012h . Here we assu me that the 
extens ion is ~0.7 of the boundaries listed in IDame et al.l 
(|1987f ) together with the mass and distances given in that 
paper. The assumption on the extension comes from the 
fact that we need to consider the 68% containment radius 
of the 7-ray emission. Wc estimate the 7-ray flux from Eq. 
[51 with k = 1, and we compare it t o the integral sen sitiv- 
ity given in (> 100 GeV, 50 hours; lActis et all 1201 ll ) and 
its scaling for the extension of the source (i.e. oc src /Sps, 
where 6>pg=0.1). The scaling of the sensitivity given here is 



The prospects for detection with CTA are slim due to 
the very large extension of the clouds. However, there is 
an abundance of smaller clouds at closer distances, that 
have been already detected at HE 7-rays. Following the 
averages der i ved fr om the cloudlet population studied in 
ITorres et al.l (|2005l) . we simulate a 40 M Q cloud at a dis- 
tance of 150 pc and an extension of 2.8 pc (1.08° in an- 
gular size). The resulting density of the cloudlet is similar 
to that of the larger cloud studied in the paragraph above, 
?ih = 140 cm~ 3 . Therefore we expect total penetration of 
CR sea for energies above E p > 1 GeV for such small radial 
extensions. However, it is unlikely to have such cloudlets 
hosting a powerful accelerator. So we consider them also as 
passive clouds. With this assumption, the expected signal is 
given by Eq. [5J with k = 1 . This translates for the cloudlet 
considered into F(> ITeV) - 2 x 10 -15 cm -2 s _1 - 0.1%o 
CU, where CU stands for Crab Units and represents the 
integral flux of the Crab nebula above 1 TeV. This flux 
would not be easily detectable (the expected CTA sensitiv- 
ity is of ~ 1%oCU). Nonetheless, a stacking approach on a 
population of roug hly 100 clouds -a list of which is used by 
ITorres et al.l (|2005f) - would lead to a detection also in VHE 
7-rays, implying a good test for the isotropy of the CR sea 
in the local neighborhood. 
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6. Peaked density profile of the molecular cloud 

In this section we describe qualitatively the expectation 
for a different density distribution of the target material 
inside the cloud. Indeed, the density of the molecular cloud 
might not be constant throughout its extension. Therefore 
we study the case of a cloud with a density profile of the 
form 



n H (r) 



n 



(1 + r/RcT 



(9) 



where Rc is the core radius and no is the densi ty at the 
center of the cloud. Following the observations bv lCrutcherl 
(fl99l we assume that the cloud magnetic field scales with 
density as 



B(r) 



100 (^lV /2 . 

V10 4 cm- 3 J 



(10) 



This will in turn affect the value of the CR diffusion co- 
efficient. We parametrize the effect of the magnetic field 
as: 



D{E,r) = D 



E/GcV V 
B(r)/3»Gj 



(11) 



The mass of the interacting target material can be calcu- 
lated following 



M 



4irm p r nn(r)dr 



47rm p noi?c 



„3-a n 



2p 2 



1 — ttr, 



(12) 



where p = (1 + r/Rc). The case of a n = reduces to a flat 
density profile investigated in Section [4] The total mass M 
is obtained with Ri = and Ri = -Rcioud- If the mass of 
the cloud is known, the central density can be derived from 
the formula above. 

We again assume M = 1O 5 M and R max =20 pc. From 
this assumption we can calculate the central density de- 
pending on the cv n chosen (with the core radius i?c=0.5 pc 
fixed). This allows us to calculate for each radius, tih, B, D 
from Eq. ©. 

With increasing densities, and correspondingly stronger 
magnetic fields, diffusion becomes much slower and the 
timescales for pp losses faster, resulting in a slower pen- 
etration of the CR sea. Because of this, the penetration of 
the CR sea into the cloud is not complete, and we can- 
not any longer assume that the sea is constant through- 
out the cloud. To compute the level of penetration of 
the CR sea into th e cloud , we follow here the approach 
by iProtheroe et alj (|2008f ) and treat the penetration as 
an analog of optical thickness, where the CR intensity 

J CR = e ~ T " ( - E ' r)l CR( E > r ') and T *( E ' r ) = T a(T a + r s ) can 
be decomposed in two terms, analogs of absorption and 
scattering: 



2Knn(r)a pp dr 



2Ka pp n Rc (1 - c*n) 1 [p 1 an ] r , 



(13) 



3D(E,r 
cR c 



■dr 



3D B Q E s p 



1 



2 



P 



(14) 



Examples of CR penetration as a function of energy are 
shown in Fig. 1101 Comparing the penetration factors and 
density profiles, we will expect a dominant contribution 
from the accelerator in the inner shells. The low number 
density of target material in the outer shell and the corre- 
spondent higher diffusion coefficient assure full penetration 
of the CR background. On the other hand, only a fraction 
of the accelerated particles will manage to diffuse out to 
the outermost shells. Indeed, because of the higher densi- 
ties found in the core of the cloud, the dominant emission in 
the outer shells will come from the CR background. With 
these assumptions, we expect to detect a source with a steep 
surface brightness profile. 
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Fig. 10. Cosmic-ray penetration factor, e~ T * t - E ' r \ for a n =2 
as a function of the particle energy for the center (R=0.5pc, 
dashed), the middle (R=10pc, solid), and the border of the 
cloud (R=19.5pc dotted), in the case of Dio = 10 26 cm 2 
s- 1 . 



7. Concluding remarks 

CTA is the forthcoming array of IACTs, and one of its 
most important physics goals is the study of the origin and 
propagation of CR. We have investigated its capabilities for 
the study of CR diffusion in molecular clouds, by making 
a phase space exploration of different cases which could be 
observable w ith such faci l ity. T his complements the study 
presented bv lAcero et al.l (|2012h . We have showed theoret- 
ical predictions for VHE 7-rays fluxes from an accelerator 
inside a massive cloud. The simulations using the CTA re- 
sponse for 50h observation time indicate that it will be pos- 
sible to constrain the diffusion coefficient parameter space. 
To have the accelerator inside the cloud is an idealized case. 
However, it allows us to study the impact of CTA observa- 
tions on expected spectral and morphological features due 
to an overdensity of CR in the cloud over the CR sea. 

We expect to be able to detect 7-ray emission from mas- 
sive molecular clouds with CTA observations. Specifically, 
for clouds of 10 5 M Q and r? src ~ 0.5°, we expect detection 
even in the case of passive cloud, i.e. only permeated by the 
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CR background, out to distances of < 1 kpc. Clouds with 
the same mass, but with an enhanced population of CR 
(k > 20) could be detected throughout the galaxy. GMC in 
the Gould Belt are very massive but very extended objects, 
therefore rendering problematic the prospect of detection 
for these clouds. 

Spectral features will aid derivation of constraints, espe- 
cially the possible identification of breaks. Moreover, emis- 
sion from extended objects could be divided in subregions. 
The superior angular resolution and sensitivity of CTA, 
with respect to the current generation of instruments, will 
permit the study of spatial spectral evolution with increas- 
ing distance from the accelerator, which retains the foot- 
print of the underlying particle distribution and its diffu- 
sion. Regarding the source profiles, older sources are ex- 
pected to have larger extension. A faster diffusion coeffi- 
cient will have the same effect. It is expected that, if the 
observed molecular cloud exhibits a peaked density profile, 
the emission profile will also be peaked. Indeed the outer 
shell of the molecular cloud might fall below detection with 
CTA sensitivity, and thus prevent the detection of the full 
extension of the cloud. Therefore the knowledge of the prop- 
erties of molecular clouds is important to aid the study of 
the objects described here. Information on the distribution 
of the target material will be of great help in order to re- 
construct the parameters relative to the acceleration and 
diffusion of charged particles in the cloud in those cases. 

At the beginning of operation, CTA might be used to 
survey the galactic plane. Two types of tests could be per- 
formed, even with small integration times for a given po- 
sition in the sky. The observation of a typical impulsive 
source like a SNR hosted in a molecular cloud with simi- 
lar mass and distance as the ones investigated here, even 
if the age of the accelerator is of the order of 10 5 years, 
will constrain the diffusion coefficient (provided a precise 
knowledge of its distance) , with detection assured for slow 
diffusion even with a few hours integration. While a source 
with the high fluxes investigated here might be rare (and 
would probably be already detected in TeV even though 
with much less detail in the spectral reconstruction), we 
have shown in Section 14.31 that we will be able to probe 
a significant fraction of the galaxy for even lower fluxes. 
Therefore, at the completion of the galactic plane survey, 
we will be able to study a numerous population, thus con- 
tributing to the study of the properties of accelerators and 
the propagation mode of CRs in molecular cloud. 

Appendix A: Grid construction 

The sphere of emission has been divided in a regular carte- 
sian grid (x,y,z). The approximation is checked with the 
cal culation of the total mass, c orrectly reproduced. Also 
the lAharonian fc Atovanl (|1996l ) curves arc reasonably re- 
produced by our code. Throught the paper the parameter 
space investigated is limited in a grid formed by the follow- 
ing parameters: 

— Diffusion coefficient: Slow to fast (e.g., Dio = 
10 26 ,I0 27 ,I0 28 cm 2 s" 1 ); 

— Diffusion coefficient energy dependence: 5 = 
0.3,0.4,0.5,0.6]; 

— Age of the accelerator: 10 3 , 10 4 , 10 5 years 

— Type of accelerator: Impulsive / Continuous 

— Spectrum of injection: 7= 2.0, 2.1, 2.2, 2.3, 2.4, 2.5; 



— Fraction of energy in input (W p = ?yl0 50 erg): t] = 
0.3,1,3. 

The corresponding flux is calculated for each grid block 
from the correspon ding proton spectrum solvin g the diffu- 
sion equation as in (EE aroman fc AtovanlfT996h . The emis- 
sion is then integrated along the line of sight. The integra- 
tion is on the variable y, where a condition is posed. The 

maximum accepted is y m ax = v -^cioud — z2 ~ x2 ■ Only a 
quarter of the sphere is calculated, because of isotropy. 

The emission integrated along the line of sight can then 
be used to simulate the flux of the chosen projected shells. 
Throughout the paper we divide the projected emission re- 
gion in five concentric shells. The projected mass of each of 
the five shells is shown in Table IA.1I 

Table A.l. Characteristics of the shells 



Rin 


Rout 


size 


MaSSp ro j 


(PC) 


(pc) 


(deg) 







4 


0.23 


6655 


4 


8 


0.23 


18590 


8 


12 


0.23 


27736 


12 


16 


0.23 


32571 


16 


20 


0.23 


23265 
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